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A new thermal model called Simple-II was presented based on modification of the original Simple 
analysis. First, the engine was modeled considering adiabatic expansion and compression spaces, in 
which effect of gas leakage from cylinder to buffer space and shuttle effect of displacer were imple¬ 
mented in the basic differential equations. Moreover, non-ideal thermal operation of the regenerator and 
the longitudinal heat conduction between heater and cooler through the regenerator wall were 
considered. Based on the magnitudes of pressure drops in heat exchangers, values of pressure in the 
expansion and compression spaces were corrected. Furthermore, based on the theory of finite speed 
thermodynamics (FST), the corresponding power loss due to the piston motion and also the mechanical 
friction were considered. Simple-II was employed for thermal simulation of a prototype Stirling engine. 
Finally, result of the new model was evaluated by comprehensive comparison of experimental results 
with those of the previous models. The output power and thermal efficiency were predicted with +20.7% 
and +7.1% errors, respectively. Also, the regenerator was demonstrated to be the main source of power 
and heat losses; nevertheless, other loss mechanisms have reasonable effects on output power and/or 
thermal efficiency of Stirling engines. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Stirling engines are one of the external heat transfer engines 
that can use various heat sources such as solar, bio-mass, convec- 
tional fuel, nuclear and so on 1-3]. Stirling engines operate in a 
closed regenerative thermodynamic cycle with two working spaces 
(compression space and expansion space) and three heat ex¬ 
changers including heater, cooler and regenerator. Compression 
takes place at a lower pressure level while heat is rejected to a cold 
heat sink through a heat exchanger called cooler. In addition, 
expansion process occurs at a higher pressure with heat addition 
from a heat source via another heat exchanger called heater. 
Working spaces of the engine are divided into expansion and 
compression cylinders with separate pistons and displacer. These 
pistons are used for displacing and distributing the working fluid 
between expansion and compression spaces. In modern Stirling 
engines, a regenerative heat exchanger called regenerator is 
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employed to recover heat from the working fluid passing between 
expansion and compression spaces. 

Developing the Stirling engine for high thermal efficiency pro¬ 
duction of power has been the concern of many researchers around 
the world. Many studies have been performed for thermal 
modeling and prediction of thermal performance of Stirling en¬ 
gines, in which researchers have tried to present precise thermal 
models that enable them to predict thermal efficiency and output 
power of Stirling engines with the highest possible precision and 
low computational cost. These researches can be classified into 
empirical, analytical and numerical methods. 

Several works can be cited here as examples for empirical 
models [4-8 . West [4] considered an output power of the Stirling 
engine and presented a new dimensionless number called Beal 
number. This number was used to predict output power of Stirling 
engines. Kongtragool and Wongwises [5] performed experimental 
analysis of Stirling engines and found that Beal number could not 
be used for predicting thermal performance of the engines working 
with high temperature ratios. Thus, they corrected Beal number 
formulation using a suitable correction factor. Prieto et al. [6-8] 
performed comprehensive experiments on Stirling engines and 
used similarity and dimensionless analysis method to create 
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empirical correlations for predicting power and efficiency of Stir¬ 
ling engines. 

Analytical models have been employed by a number of re¬ 
searchers as an option for predicting thermal performance of Stir¬ 
ling engines or Stirling cycles [9-16 . These analytical models are 
divided into two groups including methods for analyzing Stirling 
cycles [9-13] and methods for simulating Stirling engines 13-16]. 
The methods used for analyzing Stirling cycle are simple in appli¬ 
cation; but, they do not contain complicated engine parameters 
(such as dead volume, stork, heat exchanger data and so on). 
Therefore, they cannot be applied to the engine simulation in 
practical cases. These methods can be employed for thermody¬ 
namic analysis of Stirling cycles that approximately model a real 
engine. Finite time thermodynamics (FTT) and finite speed ther¬ 
modynamics (FST) are two kinds of closed form analyses that are 
employed for Stirling cycles. Yaqi et al. applied FTT for isothermal 
analysis of solar Stirling engines and tried to optimize the working 
space temperature [9]. Later, Ahmadi et al. [10] used FFT model to 
optimize a solar Stirling cycle by simultaneously considering ther¬ 
mal efficiency and output power as objective functions. Petrescu 
et al. 11 and Costa et al. 12] performed finite speed analysis of the 
isothermal Stirling engine cycle and evaluated irreversibility of the 
Stirling engine cycle in terms of speed and friction of piston 
displacement. They considered adiabatic analysis as a basic ideal 
analysis and employed experimental relations for the irreversibility 
and pressure drops in the regenerator. Also, new equations were 
used to calculate the friction factor and pressure drops in heat ex¬ 
changers. Terms of finite speed pressure drop and mechanical 
friction losses were applied as a different form in numerical model 
to calculate the power losses due to these terms as internal irre¬ 
versibility. Ahmadi et al. 13] optimized a Stirling cycle for 
maximum efficiency and power and minimum pressure loss based 
on FST model. 

The second category of analytical methods is the methods used 
for simulating real engines. These models consider engine param¬ 
eters including volume, stork, crank rotation angle, phase differ¬ 
ence between cylinders, drive mechanisms, parameters of heat 
exchangers and so on. Therefore, these models are more successful 
in simulating real engines, but are more complicated than cycle 
based models. Schmidt isothermal model 14] was the first 
analytical model introduced for determining the pressure distri¬ 
bution and performance prediction of Stirling engines. In Schmidt 
model, it was assumed that expansion and compression processes 
were isothermal. Toda et al. 15] modified Schmidt theory to 
include effect of drive mechanism type in the model and analyzed a 
gamma type Stirling engine based on the developed model. Their 
results implied that type of mechanical drives had a strong influ¬ 
ence on thermal performance of Stirling engines. Martini [16] 
presented a new mathematical model based on Schmidt model, 
in which effects of pressure drop and heat loss were considered in 
predicting output power of Stirling engines. 

Some numerical models based on thermodynamic principles 
have been considered for analysis and performance prediction of 
Stirling engines. Finkelstein 14] , for the first time, developed a new 
model based on adiabatic expansion and compression processes. In 
this model, heat exchangers were assumed to work in an 
isothermal condition. Urieli and Berchowitz [17] tried to develop 
the adiabatic model of Stirling engines using numerical calculation. 
In their work, adiabatic model was numerically calculated and a 
computer code was developed for performance prediction of Stir¬ 
ling engines. Urieli and Berchowitz [17] developed a new model 
(so-called Simple analysis) to improve results of the adiabatic 
model. In their model, non-ideal heat transfers of exchanges were 
considered. Abbas et al. [18] developed a new analysis based on 
Simple model, in which effects of the non-ideal regeneration, 


shuttle loss and heat conduction losses were considered. Effects of 
pumping work loss and regenerator heat loss were evaluated in a 
research conducted by Strauss and Dobson [19 . Granodos et al. 20] 
and Nepveu et al. 21] developed a new model based on quasi¬ 
steady flow approach for analyzing the SOLO 161 Stirling engine 
with alpha configuration. Parlak et al. [22 applied quasi-steady 
flow analysis for evaluating performance of a gamma-type Stir¬ 
ling engine. Tlili et al. [23] developed a second-order quasi-steady 
flow mathematical model for evaluating performance of medium- 
temperature differential Stirling engine. In their analysis, effects 
of shuttle loss, heat and pressure losses were considered. Effect of 
gas spring hysteresis losses in working spaces of Stirling engines 
was studied by Timoumi et al. [24 . In his analysis, effect of heat 
absorption and rejection between compression and expansion 
spaces with displacer movement (shuttle effect) was considered. 
Schulz et al. 25] presented a new general model of Stirling engines 
to predict power and efficiency. Ercan Ataer and Karabulut [26] 
used a mathematical model for the evaluation of a V-type 
Stirling-refrigerator without pressure drop analysis. Martaj et al. 
[27] analyzed a low-temperature differential Stirling engine at 
steady state operations. Robson et al. [28] presented a new math¬ 
ematical analysis of a low-temperature differential Ringbom Stir¬ 
ling engine, in which thermal parameters in expansion, 
compression and regenerator space were considered. 

Some researches have considered detailed design and evalua¬ 
tion of Stirling engines [29,30 . Cheng et al. [29] optimized the 
rhombic drive mechanism used in beta-type Stirling engines based 
on dimensionless analysis. Li et al. [30] considered effect of 
compact porous-sheets heat exchangers on performance of Stirling 
engine. In this work, heater, regenerator and cooler were formed by 
hundreds of porous metal sheets and identified for theoretical 
analysis to facilitate the future scale-up design. 

In the present work, a new numerical thermal model called 
Simple-II model was presented, which was developed by modifying 
the original Simple analysis 17] previously presented for thermal 
simulation of Stirling engines by Urieli and Berchowitz. Similar to 
the original Simple analysis, working spaces were considered as 
adiabatic spaces. The system of differential equations of the Simple 
(Adiabatic) model was modified to include effect of heat absorption 
and rejection between expansion and compression spaces by dis¬ 
placer (known as shuttle conduction heat loss) and mass leakage 
between working and buffer spaces. Therefore, differential equa¬ 
tions of mass and energy conservations of the original Simple 
analysis were modified to include terms of shuttle effect and mass 
leakage, respectively. Moreover, new forms of conservation equa¬ 
tions were used to modify the differential equation related to 
pressure. The new corrected system of differential equations was 
numerically solved using fourth-order Runge-Kutta method 17]. 
Then, similar to the original Simple model, results of the numerical 
solution were corrected for non-ideal heat transfer and pressure 
drops in heat exchangers (heater, cooler and regenerator). In final 
modification of the method, the results obtained in the previous 
parts of calculation were modified for losses due to finite speed of 
the pistons (working pressure over two pistons was corrected), 
mechanical friction and longitudinal conduction loss between 
heater and cooler through the regenerator wall. For correcting 
working pressure on two pistons, method of finite speed thermo¬ 
dynamic was applied. The developed model was applied to the 
GPU-3 Stirling engine as a case study. This engine is a prototype 
3.0 kW Stirling engine developed by General Motors. Thermal ef¬ 
ficiency and output power predicted by the Simple-II thermal 
model were compared with those of the previous thermal models 
as well as experimental data reported by Lewis Research Center 
[31]. It was shown that Simple-II model could predict thermal 
performance of Stirling engines with less deviation from actual data 
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compared to the previous thermal models. Finally, effects of engine 
parameters on thermal performance of the engine were studied 
based on Simple-II thermal model. In summary, it can be said that 
novelties of this paper were in integrating various loss effects 
including shuttle effect, mass leakage effect, power losses due to 
finite speed of pistons and mechanical friction and longitudinal 
conduction heat loss through the regenerator wall into the original 
Simple analysis, which was previously presented for thermal 
simulation of Stirling engines. 

2. Modeling 

2.2. Basic adiabatic and simple analysis 

Adiabatic expansion and compression spaces were the core of 
adiabatic and simple analyses. Adiabatic model was developed 
based on the following assumptions: 

1. All processes are performed in a steady-state condition. 

2. The engine is working at constant angular velocity. 

3. The working fluid is ideal gas. 

4. The engine is working with adiabatic compression and 
expansion spaces. 

5. Instantaneous pressure of compression and expansion 
spaces are uniform. 

6. Working fluid temperatures in the cooler and heater are 
constant. 

7. Temperature of the working fluid in the regenerator is line¬ 
arly changed. 

8. The kinetic and potential energies of gas streams is 
negligible. 

9. Heat is transferred to the working fluid only the in the heater 
and cooler. 

10. Total mass of the working gas is constant. 

11. Heat leakage between compression and expansion spaces 
and heat transfer to the environment are negligible. 

12. Leakage of gas to the outside of the engine is assumed to be 
zero. 

In the Adiabatic model presented by Urieli and Berchowitz 17], 
the Stirling engine was divided into five compartments, as illus¬ 
trated in Fig. 1. In this figure, single suffixes c, k, r, h, e represent five 
compartments including compression space, Kooler (cooler), 
regenerator, heater and expansion space, respectively. Double suf¬ 
fixes (ck, kr, rh, he) are dedicated to four interfaces between the 
compartments. 

The major governing equations were derived based on state 
equations of the ideal gas and conservation equations of mass and 
energy in each compartment. By developing the governing equa¬ 
tions for five compartments of the Stirling engine, a system of or¬ 
dinary differential equations was obtained. Table 1 indicates system 
of ordinary differential equation, mass flow equations between 
compartments and conditional temperature of the entire Stirling 
engine. In the equations of able 1, the variables including V c , V e , 


Table 1 

System of ordinary differential equations 

in the Adiabatic analysis [17]. 

p _ MR 

Vc, V k,V T V h ,V e 

T c ' T k ' T r _r T h ' T e 

Pressure 

yp ( d T Vc +fej 
£[p _ V T Ck The J 

r V/c +y(r k +r r +r h W e 

'ck \ r k ' r T h ) 'he 

Pressure 

variation 

m i = WTr 1 = C ’ k - r ’ e 

Mass 

PdV e +V e f 
dm e = fiThe ’ 

Mass 

PdV c +V c f 

dm c — RTd< 

variation 


dm,- = f, i = k,r,h 


rh ck = -dm c 
m kr = rh ck - d m k 
m rh = m he + dm h 

Mass flow 

if rh ck > 0 then T ck = T c else T ck = T k 

Conditional 

if m he > 0 then T he = T h else T he = T e 

temperature 

die = T e (f + ^-^) 

Temperature 

r\T _ t (dP,dV c dm c \ 

ai c - Ici^p + Vc mc J 

variation 

dQk - Vk R dP + Cp(T ck dm c T kr (dm c + dm k )) 

Heat output 
of cooler 

dQr = + Cp(T kr (dm c + dm k ) - r rh (dm c + dm k + dm r )) 

Heat receiving 
and rejecting 
from the 

dQn = ^^ + c p (r rh (dm c + dm k + dm r )-r he (-dm e )) 

regenerator 
Heat input 
to heater 

dW e = PdV e 

Work in 

expansion 

process 

dW c = PdV c 

Work in 

compression 

process 


dV c and dV e are obtained based on the engine configuration as 
analytic functions of the crank angle ( 0 ). \4, V Y , V h are dead volumes 
of the engine in cooler, regenerator and heater, respectively. These 
parameters are determined based on geometric specifications of 
heat exchangers and regenerator. 

Classical fourth-order Runge-Kutta method is usually used to 
solve systems of boundary-value ordinary differential equation. 

In the ideal adiabatic model, regenerator is assumed to work in 
an ideal regeneration process. In the ideal regenerative process, 
heat absorption of working fluid in the heating mode is totally 
transferred to the fluid in the cooling mode. Furthermore, pressure 
losses in the regenerator and heat exchanger are neglected. Simple 
analysis or actual Adiabatic analysis is developed to include inef¬ 
fectiveness of a real Stirling engine into the thermal model [17 . In 
Simple analysis, imperfect regeneration analysis is performed to 
determine effect of imperfect regeneration process on the magni¬ 
tude of gas temperature. Moreover, power loss due to pressure 
drops in the regenerator and heat exchangers is determined and 
subtracted from the power output obtained from the ideal Adia¬ 
batic analysis. 



Wc 


P 


Compression Kooler (Cooler) Regenerator Heater Expansion 



Fig. 1 . Schematic model for five compartments of the Stirling engine [15]. 
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2.2. Simple-II: new thermal model 

Simple-II works based on similar concepts to those of Simple 
analysis. In this model, first, effects of shuttle heat transfer loss and 
gas leakage to crank case were implemented in a differential form 
into the ideal adiabatic analysis. In this regard, the system of or¬ 
dinary differential equations of adiabatic analysis 17] was modified 
to include shuttle effect heat transfer loss and gas leakage and a 
new system of equations was obtained. In the next step, in a similar 
manner to the Simple analysis, results of the previous step were 
modified in Simple-II thermal model to include effects of pressure 
loss in the regenerator, heater and cooler, pressure loss due to finite 
motion of piston, power loss due to mechanical friction, conduction 
heat loss in the regenerator wall and non-ideal heat transfer of heat 
exchangers. It should be noted that the original Simple analysis [17] 
contained only effect of non-ideal heat transfer and pressure drop 
in heat exchangers. In the Simple-II thermal model presented in 
this paper, other loss effects (shuttle effect, gas leakage, pressure 
loss due to finite motion of piston, power loss due to the mechanical 
friction, conduction heat loss in the regenerator wall) were also 
integrated into the model for the first time. In summary, it can be 
said that, in the developed new thermal model (Simple-II), loss 
effects of Stirling engines were divided into three categories. One 
part of loss mechanisms including gas leakage and shuttle effect 
was considered in the basic ordinary differential equations. The 
second category of losses was non-ideal heat transfer and pressure 
drops in heat exchangers that were differentially evaluated in 
combination with the first part and used to correct temperature of 
engine spaces. The third category of losses was pressure loss due to 
finite motion of piston, mechanical friction and longitudinal heat 
conduction between heater and cooler through the regenerator 
wall which were calculated as separate loss terms which did not 
affect temperature distribution of engine spaces. Fig. 2 summarizes 
the three categories of loss mechanisms in an Stirling engine. 

In Simple-II thermal model, the list of assumptions of original 
Simple analysis 17] was modified and assumptions #10, #11 and 
#12 as cited in Section 2.1 of this paper for original Simple analysis 
were omitted in developing the Simple-II thermal model. 


2.2.1. Modifying ideal adiabatic analysis 

In the new model (Simple-II) in a similar manner to the ideal 
Adiabatic analysis, the engine was divided into five compartments 
of heater, Kooler (cooler), regenerator, expansion and compression 
compartments, as in Fig. 1. Also, heat transfer in the heat ex¬ 
changers was unsteady and temperature in the heater and cooler 
varied with respect to the crank angle. In developing the system of 



differential equations for governing equations, derivation of each 
parameter with respect to crank angle ( 0) was denoted by d; 
therefore, d = d/d# hereinafter. In the basic analysis of Stirling 
engines in the Simple-II model, the energy conservation law for 
each compartment of the engine (as in Fig. 1) was developed in a 
differential form to include effect of the shuttle heat loss. For each 
compartment, a hypothetical control volume which is illustrated in 
Fig. 3 was considered. For a control volume as illustrated in Fig. 3, 
equation of the first law of thermodynamics in a differential form 
was: 

dQ - dQ shuttle + (rhiCpTi - m 0 c p T 0 ) = d W + c v d(mT ) (1) 

where dQ. is heat absorption or rejection from (or into) the working 
fluid via the external resources and d W shows the net output work. 
dQshuttle denotes heat loss due to shuttle conduction by the dis¬ 
placer. This is the new term that was added to the energy conser¬ 
vation equation in the new model (Simple-II) presented in this 
paper. dQshuttle could be obtained as follows 23]: 

dQshuttle = ! ^^ L (Te-Tc) (2) 

where S, kg, Dd J, Ld, T e and T c are stork, thermal conductivity of gas, 
diameter of displacer, gap between displacer and cylinder, displacer 
length, gas temperature in expansion space and gas temperature in 
compression space, respectively. 


Heat Input to the Stilling Engine 



Fig. 2. Various loss effects of Stirling engines. 
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The third term on the left hand side of Eq. (1) determines the net 
absorbed enthalpy of the working fluid in each compartment of the 
engine due to in/out flows. The second term in the right hand of the 
energy conservation equation shows change of differential internal 
energy of the working fluid. 

Equation of state of the ideal gas for the working fluid can be 
written as: 

PV = mRT (3a) 


Assuming the adiabatic compression process, the energy con¬ 
servation equation (Eq. (1)) for compression space can be re¬ 
written as follows: 

-dQshuttie - m c k CpT C k = dWc + c v d(m c T c ) (10) 

The expression for variation of compression gas mass (dm c ) was 
obtained using state equation of ideal gases and some intermediate 
calculations as follows: 


Differentiating Eq. (3a) leads to: 

d P dV _ dm dT 

~P + ~V~ = ~m + T 

Total mass of the working fluid (M) is: 

M = m c + m k + m r + m h + m e - m leak 


dQshuttie T- CpT C kdm c — PdV c + c v d(m c Tc) 

(3b) dQ. S h Utt i e + c p T ck dm c = PdV c + c v d^-^-^j 

dQ S h U ttie + c p T ck dm c = PdV c + —^-dV c H dP 


where mi ea k denotes mass leakage to crank case due to the leakage 
effect. This term is a new factor that was added to the mass con¬ 
servation law of the original adiabatic analysis in Simple-II model. 
The leakage to the crank case was calculated as follows [17]: 


™ ^ buffer h, j J 

m leak = - ^dt- u pJ 


ARP 


g 


6/i 


P ^bufferA 

L J 


(5) 


where J, D, T g , L, u, fi, P and Pbuffer are annular gap between piston 
and cylinder, piston diameter, gas temperature, piston length, pis¬ 
ton linear velocity, viscosity of working flow, pressure of working 
space and pressure of buffer space (Pbuffer was assumed 100 kPa in 
this paper). 

With these definitions, Eq. (4) can be rewritten in the following 
form: 


M = 


pfYc_ k v k v L ,v k ,Ve 
1 [Tc + T k T r + T h Te 


- < 7rD 


R 

P "E ^buffer 


(“*> - fe 


( 6 ) 


4RT g 6/d 


The mass equation in a differential form was obtained by 
differentiating Eq. (4) as follows: 


d in c + dm k + d in r + dm h + dm e - dm leak = 0 


(7) 


As for the heater, cooler and regenerator, the following was 
assumed: dV/V = 0. 

Therefore, differential mass equation (Eq. (3b)) can be written as 
follows: 


dP 


dm 


k,h,r 


Therefore, 


m 


k,h,r 


(8a) 


dm 


mdP 


k,h.r 


VdP 


k,h,r 


RT 


k,h,r 


(8b) 


Substitution of Eq. (8b) for Eq. (7) led to the following expression: 


dQshuuie + c p T ck dm c = ^PdVc+^VcdP 


dQ< 


+ r c|( dm 


PdVr C\ 

+ 


R Rc 


VrdP 


DQshuttie 

Qk C P 


+ dm c 


^ + _L %p 

RT ck RT ck 7 


PdVc + v fdP d Q shuttle 


d 777 r 


RT, 


ck 


c p Qk 


( 11 ) 


In a similar manner for the expansion space, the following can 
be obtained: 


dmp = Pdl4+!dP_dQ^ 

Cp 1 he 


PT, 


( 12 ) 


he 


Substituting dm c and dm e obtained from Eqs. (11) and (12) for 
Eq. (9) led to the following expression of differential pressure in the 
working space: 


dP = 


-P[ dVc + + dQshuttleft Mhe + Pdm leak 


ck 


' he 


14 | (Yk I Vh | VA I V e 
jT ck ^ \T k + T h + rj + yThe 


(13) 


Moreover, the differential temperature equations were obtained 
as follows (see Eq. (3b)): 


, T f dP dV c 

dTc = 7 % + V7 

dP dVe 


dip = Te 


+ 


V e 


dm c 

m c 

dm ( 

m P 


(14) 


Heat transfer rate in the cooler, heater and regenerator could be 
evaluated from the energy conservation equation as follows: 


dQk 


V k dPc v 

R 


Cp ( T ck m ck T kr m kr ) 


dQh = - c P (T he m he - T he m he ) 


(15) 


(16) 


dm c + dm e - dm leak + 



R 


0 



dQr 


V r d Pc v 
R 


c p (T kr m kr - T rh m rh ) 


(17) 
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The output work performed by the power piston can be calcu- 


lated as follows: 


dW c = PdV c 

(18) 

dW e = PdV e 

(19) 


2.2.2. Implication of other loss effects considered in Simple-II model 
System of ordinary differential equations obtained in Section 
2.2.1 was solved using the fourth-order Ruge-Kutta method in a 
similar manner that differential equations were solved in original 
ideal adiabatic analysis 17]. The difference between the present 
model and original ideal adiabatic model was that, in the present 
model, shuttle effect and gas leakage terms were included in 
differential equations. Effect of the losses including non-ideal heat 
transfer in heat exchangers, pressure drops in heat exchanger, 
pressure loss due to finite motion of piston, power loss due to 
mechanical friction and conduction heat loss in the regenerator 
wall was divided into two groups. The first group including non¬ 
ideal heat transfer and pressure drops in heat exchangers was 
differentially evaluated in combination with the first part and 
used to correct temperature of the engine spaces. The next group 
of losses including pressure loss due to finite motion of the piston, 
mechanical friction and longitudinal heat conduction between the 
heater and cooler through the regenerator wall was calculated as 
separate loss terms that did not affect temperature distribution of 
engine spaces; these loss effects were subtracted from the 
results obtained from solution of differential equation. In the 
following sections, principles or evaluation of these losses are 
presented. 


£ 


1 


1 + 


5T 

Thi-Tho 


( 22 ) 


where 5T is temperature difference between the hot and cold 
streams in the regenerator and Thi and Th 0 are hot stream tem¬ 
peratures at the input and output of the regenerator. 

The energy balance equation for the hot stream 
( mc p (T hi - T h0 ) = hA wg 5T) leads to the following expression for 
regenerator effectiveness: 


1 


1 


NTU 


£ = 


1 | mc P 


hA 


1 + 


wg 


NTU 


1 + NTU 


(23) 


Convection heat transfer coefficient (h) can be calculated as 
follows [32]: 


St = 0.023 Re-° 2 Pr-°- 6 


(24) 


NTU 


StA 


Wg 


A 


(25) 


In Eqs. (23)—(25), NTU, St, A wg , A, h, Re and Pr are the number of 
transfer unit, Stanton number, wet area of gas in the heat 
exchanger, cross-section area of heat exchanger, convective heat 
transfer coefficient, Reynolds number and Prandtl number, 
respectively. In the evaluation of the Reynolds number, hydraulic 
diameter of the regenerator mesh is obtained as [32]: 



4U 

T>(i - n) 


(26) 


2.2.2A. Non-ideal heat transfer effect. In Stirling engines, heat is 
transferred to the working fluid in the heater and cooler. Then, it is 
absorbed and recovered in a regenerative heat exchanger, known as 
the regenerator. For perfect heat recovery, the regenerator is 
designed to be very compact and efficient. The main heat transfer 
mechanism in the regenerator is forced convection heat transfer. 
Special structure of the regenerator increases frictional pressure 
loss (pumping loss). 

Performance of Stirling engines for heat recovery is evaluated by 
regenerator effectiveness. Effectiveness is defined as a ratio of real 
enthalpy change to maximum enthalpy change of working fluid in 
the regenerator. Regenerator effectiveness changes from 0.0 (for no 
heat recovery) to 1.0 (for complete heat recovery). When effec¬ 
tiveness is less than unity, the output gas temperature from the 
regenerator is less than temperature of the heater. This heat should 
be supplied by an external heat source, which causes efficiency 
reduction. These definitions can be written for the heater and 
cooler, respectively, as follows: 


Qh — Qh ldea l + Qrioss — Qhideal + QridealO ~ £ ) (20) 


Ok — Qkjdeal “ Qrioss ~ Qkideal _ Ordeal 0 ~~ £ ) 


( 21 ) 


where Qr loss is heat loss in the regenerator, Qh is heat addition to the 
working fluid in the heater, Qk is heat rejection to the cooler and £ is 
the regenerator effectiveness. 

Regenerator effectiveness could be evaluated as a function of the 
input and output temperatures by assuming a linear profile in the 
regenerator 17]: 


<|> = Av/ s 


V, 


mesh 


(27) 


where V mes h* n and <f> are volume of wires, porosity and shape 
factor of the regenerator mesh, respectively. 

Because of non-ideal heater and cooler, the working fluid tem¬ 
perature in these two heat exchangers are lower and higher than 
the wall temperature, respectively. By calculating the heat loss in 
non-ideal regenerator and heat transfer calculation, the wall tem¬ 
perature in cooler and heater can be calculated as follows: 


n n n “hA/vh(^wh ^h) rno\ 

Oil = Qh lde al +( 2rioss = - J - US) 

where f T w h, /ih» Avh and Qh are engine rotation frequency, wall 
temperature in the heater, convective heat transfer of the working 
fluid in the heater, heat transfer area of the heater and rate of heat 
transfer in the heater, respectively. Therefore, the heater temper¬ 
ature is corrected from the value of heater wall temperature as 
follows: 




QJ 

AiA/vh 


(29) 


In a similar manner, the cooler temperature is corrected as: 




QJ 

A<A/vk 


(30) 


In Eqs. (29) and (30), heat transfer coefficient is obtained using 
the following correlation 17] : 
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0.0791 ii-Cp-Re^ 75 
2D h Pr 


(31) 


where Re m and D H are Reynolds number and hydraulic diameter of 
heater or cooler. 

Corrected values of heater and cooler temperatures give feed¬ 
back to the solver of the system of differential equations (cited in 
Section 2.2.1) until the convergence criteria for the magnitude of 
gas temperatures is satisfied. 


2.22.2. Pressure drop effect. Pressure drops in heat exchangers due 
to the fluid friction cause power loss in Stirling engines. The pres¬ 
sure drop in heat exchangers should be calculated in order to 
calculate the power loss due to pressure losses. The pressure drops 
in heat exchangers can be obtained as: 


A P 


d 2 A 


(32) 


where d is hydraulic diameter, u is velocity, V is volume, A is flow 
cross-section area and fR e is a Reynolds friction factor .f Re is defined 
as follows [17]: 



f 16 Re < 2000 

f Re = \ 7.343 x 10- 4 tfe 13142 2000 < Re < 4000 (33) 

{ 0.0791 Re 0 75 Re > 4000 

The net power loss due to pressure drop in heat exchangers of 
the Stirling engine is evaluated as follows: 

W loss = J J>Pd\4 (34) 


2.2.23. Effect of finite speed of piston and mechanical friction. 
Based on finite speed thermodynamic principle, the pressure in 
compression and expansion pistons differs from the calculated 
instantaneous pressure of expansion and compression spaces. It 
was demonstrated that pressure over the piston in expansion space 
was less than the mean pressure of the cylinder in the expansion 
process. Similarly, the pressure over the piston surface in 
compression process is greater than mean pressure of the cylinder. 
Therefore, the work obtained in expansion process is less than that 
of theoretical expansion. Similarly, the actual compression work is 
greater than the theoretical work of compression which is calcu¬ 
lated in classical thermodynamics. Therefore, output power of the 
real engine decreases from the corresponding value calculated 
using classical thermodynamics. The pressure loss caused by finite 
speed effect and mechanical friction can be calculated using Eqs. 
(35) and (36) as follows [33]: 


cylinder-corrected 


cylinder 



d= 


AP f \ 


cylinder 


) 


(35) 


AD ClUp APf / 

1X1 finite speed = ±—~ ± n- (3bJ 

L 1 cylinder 

Finite speed movement of the piston that generates pressure 
waves in working spaces causes finite speed pressure losses. In the 
above Eqs. (35) and (36), sign (+) is used for compression space and 
sign (-) are used for expansion space. Also, u p is the piston velocity, 
c is average speed of molecules and APf is piston-cylinder friction 
loss. In Eqs. (35) or (36), APf, a and c are calculated as follows [34]: 


Fig. 4. Flow-chart of the Simple-II thermal model. 


Table 2 

System of ordinary differential equations obtained in the new (Simple-II) thermal 
model. 
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Table 3 

Design specifications of the GPU-3 Stirling engine [26]. 


General information Working fluid 

Helium 

Regenerator 

Regenerator internal diameter 

22.6 mm 

Piston stroke 

31.2 mm 


Regenerator length 

22.6 mm 

Internal diameter of cylinder 

69.9 mm 


Number of regenerator 

8 

Frequency 

41.7 


Material 

Stainless steel 

Mean pressure 

4.13 MPa 


Diameter of regenerator tube 

0.04 mm 

Phase angle 

90 


Porosity 

0.697 

Heater temperature 

977 K 

Cooler 

Average of tube length 

46.1 mm 

Cooler temperature 

288 K 




Number of Cylinder 

1 


Internal diameter of tube 

1.09 mm 

Heater Average tube length 

245.3 mm 


External diameter of tube 

1.59 mm 

Tube external diameter 

4.83 mm 


Number of tubes 

312 

Number of tubes 

40 

Volume 

Swept Vol. (expansion/compression) 

120.88/113.14 mm 3 

Tube internal diameter 

3.02 mm 




Tube length (cylinder side) 

116.4 mm 


Clearance Vol. (expansion/compression) 

30.52/28.68 mm 3 

Tube length (regenerator side) 

128.9 mm 


Dead Vol. (heater/cooler/regenerator) 

70.88/13.8/50.55 mm 3 


AP f = 0.97 + 0.15^ (37a) 

a = +3? (37b) 

c = V3RT (37c) 

where N r is angular rotational speed of the engine in rpm. 

Finally, the work loss due to finite speed of the piston motion 
and mechanical friction can be calculated as follows: 

W,oss-finite speed = / (± ^ ± (38) 

J \ c 1 cylinder ) 


222.4. Effect of heat loss from the regenerator wall. Heat conduc¬ 
tion leakage between the heater and cooler through the regener¬ 
ator wall is evaluated using Eq. (39) [35]: 

Qreg-leak = ^cond(^wh — ^wk) (39) 

where R c on d» T w h and T W k are conduction resistance of the regen¬ 
erator wall and wall temperatures of the heater and cooler, 
respectively. 

3. Solution method 

Solution algorithm of Simple-II thermal model is illustrated in 

Fig. 4. 

As is clear in Fig. 4, in the first step of the calculation, system of 
ordinary differential equation obtained in Section 2.2.1 was 
numerically solved. This step was quite similar to the original ideal 
adiabatic analysis, except that system of differential equations was 
corrected to include effect of shuttle heat transfer loss and mass 


leakage. In a similar manner to able 1, system of ordinary differ¬ 
ential equations is summarized in fable 2. 

In the system of equations, variables including V 0 V e , dV c and 
dV e are obtained based on the engine configuration as analytic 
functions of the crank angle 0. \4, V Y and V h are dead volumes of the 
engine in cooler, regenerator and heater, determined based on 
geometrical specifications of heat exchangers and regenerator. 
Regenerator temperature is considered as mean effective temper¬ 
ature of heater and cooler [17 . Except for the aforementioned 
constants, it is required to solve 16 ordinary differential equations 
including 24 variables over a complete cycle (see Table 2), among 
which seven derivatives (T c , T e Qi<, Qr» Qh, W c and W e ) are numerically 
integrated using the fourth-order Runge-Kutta method [17] and 
eleven variables and derivatives including W, p, Qshuttie, V Ct V e , m c , 
mk, m r , m^, m e and mi ea k are analytically determined. Moreover, 
there are six boundary conditions and mass flow rate equations the 
same as those indicated in Table 1 for original ideal adiabatic 
analysis. In this case, the equations were initial-value problems. 
Initial estimation of the working fluid mass was performed using 
Schmidt model 14]. Initial gas temperatures in compression and 
expansion cylinders were estimated to be equal to the cooler and 
heater temperatures, respectively. Then, numerical solutions of 
differential equations were performed for several times until the 
steady-state condition was obtained, in which the magnitudes of 
temperatures at the end of cycle converged to the corresponding 
temperatures at the beginning of the cycle. In other words, 
magnitude of temperatures at the start of the cycle ((9 = 0) had to be 
equal to the corresponding values at the end of the cycle (6 = it). 
Hence, the calculations were repeated until the same gas temper¬ 
atures were obtained for 0 = 0 and 0 = 7T. 

As is clear in Fig. 4, in the next step, calculation was entered into 
the Simple analysis section, in which the non-ideal heat transfer 
effect was taken into account to correct magnitude of cooler and 
heater temperatures using Eqs. (29) and (30). The corrected mag¬ 
nitudes of cooler and heater temperatures gave feedback to the 
numerical solution section until convergence for the magnitude of 


Table 4 

Comparison between the current and other models. 


Model 

Ideal adiabatic model 

Urieli & Berchowitz Model 

Timoumi model [19] 

Current study (Simple-II 

Experimental 




(Simple Analysis) [12] 



model) 


result [31 ] 


Actual data 

Error a (%) 

Actual data 

Error 3 (%) 

Actual data 

Error 3 (%) 

Actual data 

Error 3 (%) 


Efficiency (%) 

62.3 

+40.0 

52.5 

+31.2 

38.5 

+17.2 

28.4 

+7.1 

21.3 

Power (kW) 

8.30 

+176.7 

6.70 

+123.3 

4.27 b 

+42.3 

3.62 

+20.7 

3.0 


a Error in thermal efficiency was calculated as a difference between the actual value and experimental data. 
b Hysteresis loss was considered (without the effect of hysteresis power was obtained as 5.90 kW in Ref. [19]). 
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Comparison between the current results and experimental results at various engine operation modes. 
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Frequency (Hz) Mean effective Thermal efficiency predicted by the simple analysis Thermal efficiency predicted by the current study Experimental 

pressure (MPa) (variable Pr) (Simple-II model) (variable Pr) efficiency [26] 




Actual 
value (%) 

Error 

(as differences) (%) 

Average error 
(as difference) (%) 

Actual 
value (%) 

Error 

(as differences) (%) 

Average error 
(as difference) (%) 


16.67 

2.76 

38.72 

18.22 

17.90 

32.48 

11.98 

12.85 

20.50 

25.00 

2.76 

36.16 

15.46 


31.21 

10.51 


20.70 

33.33 

2.76 

33.79 

15.79 


29.45 

11.45 


18.00 

41.67 

2.76 

31.48 

16.28 


27.45 

12.25 


15.20 

50.00 

2.76 

29.12 

17.32 


25.21 

13.41 


11.80 

58.33 

2.76 

29.74 

24.34 


22.89 

17.49 


5.40 

25.00 

4.14 

35.65 

10.85 

11.46 

32.29 

7.49 

8.28 

24.80 

33.33 

4.14 

33.52 

9.62 


30.40 

6.50 


23.90 

41.67 

4.14 

31.48 

10.18 


28.39 

7.09 


21.30 

50.00 

4.14 

29.45 

11.25 


26.33 

8.13 


18.20 

58.33 

4.14 

27.40 

15.40 


24.21 

12.21 


12.00 

41.67 

5.52 

31.20 

8.70 

10.82 

28.59 

6.09 

8.11 

22.50 

50.00 

5.52 

29.33 

10.53 


26.62 

7.82 


18.80 

58.33 

5.52 

27.44 

13.24 


24.62 

10.42 


14.20 

50.00 

6.90 

29.07 

10.37 

11.73 

26.61 

7.91 

9.19 

18.70 

58.33 

6.90 

27.29 

13.09 


24.67 

10.47 


14.20 



Frequency (Hz) 


(a) 



(b) 


Fig. 5. Comparison of thermal efficiency between original Simple analysis and the current study’s (Simple-II thermal model) at (a) various frequencies and P mea n = 2.76 MPa; (b) 
various mean effective pressures and fr = 50 Hz. 
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temperatures was obtained. Finally, magnitude of output power 
was corrected for power losses due to pressure drops in heat ex¬ 
changers, mechanical friction and finite motion of the piston. In 
addition, the magnitude of heat transfer in the heater was corrected 
for including effect of conductive heat loss through the regenerator 
wall. 

4. Case study 

Evaluation of the developed thermal model was performed by 
considering the GPU-3 Stirling engine as a case study. The GPU-3 
Stirling engine was built in research laboratories of General Mo¬ 
tors as a 3.0 kW Stirling engine and was tested in the NASA Lewis 
Research Center. This engine is a beta-type engine with the rhombic 
drive and compact heat exchangers and regenerator with a woven 
mesh. Design specifications of the GPU-3 Stirling engine are indi¬ 
cated in Table 3. 


5. Results and discussion 

5.2. Model verification and error analysis 

Comparison between the results of the current new thermal 
model (Simple-II) and corresponding results obtained from previ¬ 
ous models and experimental data were indicated in Table 4. 
Comparisons of the current results with experimental results 
showed that Simple-II model had better agreement with the 
experimental results in comparison to the previous thermal models 
so that errors in the prediction of the output power and thermal 
efficiency were reduced from 42.3% and 17.2% (as difference) in the 
best previous model 19] to 20.7% and 7.1% (as difference) in the 
present model, respectively. Therefore, the current thermal model 
improved errors on the output power and thermal efficiency by 
about 21.6% and 10.1%, respectively, in comparison to the best of 
previous thermal models [19 . 
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Fig. 6. Comparison of thermal efficiency error between original Simple analysis and the current study’s (Simple-II thermal model) at (a) various frequencies and P mea n = 2.76 MPa; 
(b) various mean effective pressures and fr = 50 Hz. 
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Similarly, ^able 5 compares Simple-II thermal model and Simple 
model in terms of predicting thermal efficiency at various engine 
rotation speeds (rotation frequencies) and mean effective pres¬ 
sures. This table indicates that, in general, error in prediction of 
thermal efficiency for both models decreased as mean effective 
pressure of the engine increased. Similarly, it can be said that the 
error of both models was generally higher when rotation frequency 
of the engine increased or rotation speed was very low. fable 5 re¬ 
confirms that the presented thermal model predicted thermal ef¬ 
ficiency with lower error in comparison to the original Simple 
analysis. Based on Table 5, thermal efficiency of the Stirling engine 
reduced with increasing rotation frequency or mean effective 
pressure. Based on Table 5, thermal efficiency of the Stirling engine 
reduced with increasing rotation frequency or mean effective 
pressure. 

For a better insight into effect on rotation frequency and mean 
pressure on thermal efficiency of the engine, Fig. 5a and b is pre¬ 
sented here, respectively. Effect of frequency on efficiency of Stir¬ 
ling engine at constant pressure (2.76 MPa) is illustrated in Fig. 5a, 
which confirms that thermal efficiency of the engine reduced when 
rotation frequency increased. It could be found in Fig. 5a that both 
models had reasonable deviation from the experimental data; 
nevertheless, the proposed Simple-II thermal model not only had 
less deviation but also showed a better matching variation trend to 
the actual experimental trend. Based on the original Simple anal¬ 
ysis, thermal efficiency decreased linearly when rotation frequency 
decreased to 50 Hz and increased linear again as rotation frequency 
increased to above 50 Hz. However, as is clear in Fig. 5a, experi¬ 
mental data demonstrated that thermal efficiency of the engine 
decreased for the entire increasing range of rotation speed. In this 
regard, the variation trend of thermal efficiency with rotation speed 
in Simple-II thermal model was quite similar to that of the exper¬ 
imental data. 

Effect of effective mean pressure on thermal efficiency of Stirling 
engine at constant rotation frequency (fr = 50 Hz) is depicted in 
Fig. 5b. It is clear that Simple-II model predicted thermal efficiency 
with less deviation from real data in comparison to the original 
Simple model. Based on the original Simple model, thermal effi¬ 
ciency remained almost constant when the mean pressure was 
increased to 5.5 MPa and linearly increased with a higher rate when 
mean pressure increased to above 5.5 MPa. Based on the present 
model, a small linear increasing rate of thermal efficiency was 


observed when mean pressure was increased over the entire vari¬ 
ation range of the mean pressure. Actual data showed that, over 
4.0 MPa, increasing rate of thermal efficiency was similar to what 
predicted by Simple-II model; but, when mean pressure was less 
than 4.0 MPa, increasing rate of thermal efficiency of the real data 
was higher than what was predicted by both models. 

Fig. 6a and b illustrates variation in errors of two models in 
terms of predicting thermal efficiency against variation of rotation 
frequency and mean pressure, respectively. Based on Fig. 6a, errors 
in the prediction of thermal efficiency for both models increased 
linearly when frequency increased from 25 Hz to 50 Hz; rate of 
increase in error increased drastically when frequency increased to 
above 50 Hz. The linear increasing rate of error for Simple-II model 
was a little higher than that of Simple analysis when frequency was 
smaller than 50 Hz. When frequency was over 50 Hz, linear 
increasing rate of Simple-II model was lower than the corre¬ 
sponding error predicted by the Simple analysis. As is clear, abso¬ 
lute value of error obtained in Simple-II model was lower than the 
corresponding error of Simple analysis for all the values related to 
engine frequency. Fig. 6b indicates that, for both models, magni¬ 
tude of error in prediction of thermal efficiency (as difference) 
linearly decreased when mean pressure increased up to =4.2 MPa. 
Over 4.2 MPa, variation errors of both models were non-significant 
as mean pressure varied. Variation rate of errors for both models 
was almost the same; however, in all case, Simple-II model had less 
error than Simple analysis. 

5.2. Performance prediction of the case studied Stirling engine 

P-V diagram of the compression and expansion spaces pre¬ 
dicted by Simple-II model is illustrated in Fig. 7. 

Pressure drops in heater, cooler and regenerator at various crank 
angles are shown in Fig. 8a. As pressure drops of heat exchangers 
were calculated in both Simple and the Simple-II models, no dif¬ 
ference in pressure drop prediction of the two models was 
observed in Fig. 8a. The results showed that the fluid friction in the 
regenerator had a significant contribution to pressure drop and 
power loss of the Stirling engine. As magnitudes of pressure drops 
in heater and cooler were much smaller than the corresponding 
pressure drop in the regenerator, on the scale of Fig. 8a, variation of 
these two pressure drops against crank angle were not clear; 
therefore, Fig. 8b is presented here to clearly show variation of 



Fig. 7. Pressure distribution at various volumes of compression and expansion spaces at fr = 41.67 Hz and P mea n 


= 4.14 MPa. 
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(a) 



Crank Angle (deg) 

(b) 


Fig. 8. (a) Pressure drops in three heat exchangers (regenerator, heater and cooler) at various crank angles at fr = 41.67 Hz and P mea n = 4.14 MPa; (b) pressure drops in heater and 
cooler at various crank angles at fr = 41.67 Hz and P mean = 4.14 MPa. 


pressure drops in the heater and cooler against crank angle. This 
figure indicates that the pressure drop in the heater was about 50% 
higher than the corresponding pressure drop of the cooler. Com¬ 
parison of Fig. 8a and b indicated that maximum pressure drop in 
the regenerator was about 8.0 times of the sum of maximum 
pressure drops of the heater and cooler. 

Pressure drop due to finite motion of pistons at various crank 
angles for compression and expansion spaces is depicted in Fig. 9. 
The results indicated that the pressure drop related to finite speed 
of the piston in compression space was higher than the corre¬ 
sponding pressure drop of the expansion space due to the fact that 
the working fluid temperature in the compression space was lower 
than the one in expansion space and finite speed pressure drop was 
in proportion to the inverse of the fluid temperature (1/T). 

Comparison of Figs. 8(a, b) and 9 revealed that, although the 
highest pressure drop occurred in the regenerator, pressure drop 
due to finite speed of the piston had a reasonable value. Maximum 
pressure loss due to finite speed of the pistons was about 30% of the 
highest pressure drop of the regenerator. Moreover, as can be seen, 


pressure drops due to finite speed of the piston were much higher 
than those of heater and cooler. For a better insight into effect of 
each loss mechanism on total loss of the GPU-3 Stirling engine, 
Table 6 is presented here. 

Based on Table 6, almost 20.6% of the generated power was lost 
due to various power loss mechanisms, in which the most part of 
the power loss was due to pressure drop in heat exchangers so that 
76.8% of total power loss (equal to 15.8% of generating power) was 
due to pressure drop in heat exchangers. The second and third 
portions of power loss belonged to finite speed of the piston and 
mechanical friction, respectively. Power loss due to gas leakage to 
the buffer space was 3.0% of total power loss as about 0.6% of the 
generated power. Moreover, as the shuttle effect was coupled to the 
developed system of differential equations, it was possible to 
evaluate its effect not only on heat loss but also for the power loss. 
Therefore, it was found that the shuttle effect was responsible for 
0.6% of total power loss equal to 0.1% of the generated power. 
However, the shuttle effect had a minor role in the power loss of 
Stirling engine (based on Table 6); nevertheless, as is clear in 
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Fig. 9. Finite speed pressure drop at various crank angles for compression and expansion spaces at fr = 41.67 Hz and P me an = 4.14 MPa. 


Table 6 

Surveying the contribution of various loss mechanisms of GPU-3 Stirling engine. 


Loss mechanism 


Percent of the total generated power or total heat addition 

Percent of the total loss 

Power losses 

Power loss due to finite speed of the piston 

3.2 

15.2 


Power loss due to mechanical friction 

0.9 

4.4 


Power loss due to pressure drops in heat exchangers 

15.8 

76.8 


Power loss due to mass leakage 

0.6 

3.0 


Power loss due to shuttle effect 

0.1 

0.6 


Total power loss 

20.6 

100.0 

Heat losses 

Heat loss due to regenerator ineffectiveness 

29.5 

78.8 


Longitudinal conduction heat loss 

4.5 

12.0 


Heat loss due to the shuttle effect 

3.0 

8.0 


Heat loss due to mass leakage 

0.5 

1.2 


Total heat loss 

37.5 

100 


Table 6, the shuttle loss was responsible for 8.0% of total heat loss of 
the GPU-3 Stirling engine (3.0% of total heat addition to engine from 
heat source). Moreover, it can be seen that, among the heat loss 
mechanisms, the greatest heat loss was due to regenerator 


ineffectiveness (78.7%) followed by the longitudinal conduction 
loss through the regenerator wall (12.0%) and shuttle effect (8.0%). 
Contribution of mass leakage to total heat loss of the GPU-3 Stirling 
engine was 1.2%. Based on able 6, it is clear that the regenerator 



Fig. 10. Temperature distributions in heater, cooler and regenerator at various crank angles at fr = 41.67 Hz and P mea n = 4.14 MPa. 
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(a) 



(b) 

Fig. 11. Effect of mass leakage to the buffer space on (a) thermal efficiency; (b) output power (fr = 41.67 Hz and P me an = 4.14 MPa). 


was responsible for greater portion of power loss and heat loss. 
Therefore, the regenerator is a bottleneck for any thermal 
improvement task of Stirling engines. The role of mass leakage to 
buffer space was considerable in power loss, but minor in heat loss. 
On the other hand, shuttle effect had significant contribution to 
heat loss, but non-significant effect on power loss. Therefore, all the 
loss mechanisms considered in this paper had more or less effect on 
thermal performance of Stirling engines. 

The temperature distribution in heater, cooler and regenerator 
at various crank angles is depicted based on Simple-II model and 
Simple analysis in Fig. 10. It can be observed that temperature of the 
expansion heater that was predicted by the Simple-II model was 
slightly less than what was predicted by the Simple analysis. On the 
other hand, the Simple-II model predicted slightly higher temper¬ 
ature in the compression space compared to the predicted 
compression temperature by the Simple analysis. Difference in 
temperature prediction by of the Simple-II thermal model and 
Simple analysis was due to consideration of effects of shuttle loss 
and longitudinal conduction loss through the regenerator wall in 
the Simple-II model. 


Effect of cylinder-piston gap in thermal efficiency and output 
power of the Stirling engine is depicted in Fig. 11a and b, respec¬ 
tively. Results showed that thermal efficiency and power decreased 
with increase in the cylinder-piston gap due to increasing mass 
leakage. At very low values of gap dimensionless number (J/D), 
decreasing rate of thermal efficiency and power against increasing 
J/D was non-significant; but, when J/D increased to above 
0.6 x 1CT 3 , thermal efficiency decreased drastically. At J/ 
D = 0.4 x 1CT 3 as design J/D of the GPU-3 Stirling engine, which is 
shown Fig. 11a and b, effect of gas leakage to the crank case on 
thermal performance of the engine was non-significant. 

Effect of the heater-cooler temperature difference on thermal 
efficiency and output power are is considered in Fig. 12a and b, 
respectively. These figures indicate that the thermal efficiency and 
output power of the engine increased with increasing the heater 
and cooler temperature difference. Moreover, increasing rates of 
thermal efficiency and power were the same by both Simple-II 
model and Simple analysis. 

As previously cited in the interpretation of Fig. 8a, the regen¬ 
erator is responsible for a great part of the system power losses due 
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(b) 


Fig. 12. Comparing effect of heater-cooler temperature difference in (a) thermal efficiency; (b) output power obtained by in Simple and Simple-II models (fr = 41.67 and 
Anean = 4.14 MPa). 


to very high-pressure drop. On the other hand, the regenerator is an 
essential part of Stirling engines as it recovers a great part of 
thermal energy of the working fluid in order to increase thermal 
efficiency of the engine. In addition, porosity of the regenerator 
mesh in a key parameter in thermohydraulic performance of the 
regenerator because pressure drop and heat regeneration of the 
regenerator are highly dependent to the porosity factor. Figs. 13a 
and b are presented here to show effect of porosity variation in 
thermal performance and output power of the GPU-3 Stirling en¬ 
gine, respectively. Based on Fig. 13a, increasing porosity of the 
regenerator (using a regenerator matrix with lower density) up to 
0.80 led to increasing thermal efficiency since a reduction in 
porosity caused reduction in power loss. On the other hand, 
decreasing the regenerator porosity caused reduction of heat re¬ 
covery. Therefore, as is clear in Fig. 13b, thermal efficiency-porosity 
curve showed a peak at <p = 0.80. Fig. 13b demonstrates that the 
output power increased when regenerator porosity decreased, 
which was due to reduction in the pressure loss in the regenerator 


at lower values of porosity. Design value of GPU-3 porosity was 
0.696. Therefore, if porosity of the regenerator increased from 0.696 
to 0.80 (the optimum value), the predicted thermal efficiency 
increased from 28.4% to 32.5% and output power increased from 
5.90 kW to 6.0 kW (based on output of simple-II thermal model). 

Fig. 14 shows variation of the power loss against rotation fre¬ 
quency of the GPU-3 Stirling engine. Both Simple-II thermal model 
and Simple analysis predicted that power loss increased rapidly 
when rotation frequency of the engine increase. Based on Fig. 14, 
two graphs diverged when rotation frequency increased, which 
was due to the fact that, at higher values of engine frequency, po¬ 
wer loss increased due to finite speed of the piston and mechanical 
friction. Therefore, deviation increased between power losses 
predicted by two models at higher frequency. At lower values of 
engine frequency, as power losses due to finite speed of the piston 
and mechanical friction became non-significant, the power losses 
predicted by Simple-II thermal model and Simple analysis got close 
to each other. 
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(a) 



Matrix Porosity 

(b) 


Fig. 13. Effect of regenerator matrix porosity on (a) thermal efficiency; (b) output power obtained at fr = 41.67 Hz and P mea n = 4.14 MPa. 


6. Conclusions 

The famous numerical thermal model of Stirling engines, called 
Simple model, was modified and a new thermal model called 
Simple-II, which considered various loss mechanisms including 
power loss due to finite speed of the piston, mechanical friction, 
shuttle effect, longitudinal conduction loss between heater and 
cooler through the regenerator wall and mass leakage, was devel¬ 
oped. Based on the methodology presented in this work, the system 
of ordinary differential equation obtained in the analysis was 
solved using fourth-order Runge-Kutta method and results of the 
numerical solution were corrected for effects of various losses. 
Simple-II model was employed for thermal simulation of a proto¬ 
type Stirling engine called GPU-3 engine. It was shown that Simple- 
II thermal model yielded results that were more consistent with 
those of the experimental data in comparison to the previous 
thermal models. It was also demonstrated that variation trends of 
thermal efficiency against variation of the engine rotation fre¬ 
quency and mean pressure of the Simple-II thermal model were 


more consistent with the experimental trends compared to the 
original Simple analysis. It was observed that the prediction error of 
thermal efficiency was higher at higher engine rotation speeds and 
lower mean pressures. Based on the comprehensive analysis, it was 
concluded that regenerator was responsible for the greater part of 
power and heat losses due to pressure drop and thermal ineffec¬ 
tiveness, respectively. After the regenerator, the following parts or 
effects were responsible for the greatest values of power loss in a 
descending order: finite speed of piston, mechanical friction, heater 
pressure drop, cooler pressure drop, gas leakage and shuttle effect. 
Moreover, after the regenerator, the greatest part of the heat loss 
was due to longitudinal conduction loss through the regenerator 
wall followed by shuttle effect and mass leakage, respectively. It 
was concluded that, if the piston-cylinder gap is in a reasonable 
range of gaps used for designing engines, the leakage effect would 
have a minor effect on performance of the engine; but, at large 
gaps, the output power and thermal efficiency would drastically 
reduce. Porosity analysis of regenerator indicated that, if the cur¬ 
rent porosity of the regenerator mesh of the GPU-3 Stirling engine 
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Fig. 14. Effect of rotation frequency on power loss at P me an = 2.76 MPa. 


increased by about 14.8%, the corresponding thermal efficiency and 
output power of the engine would increase by 4.1% (as difference) 
and 1.7%, respectively. As effect of power loss due to finite speed of 
the piston and mechanical friction was more significant at higher 
engine rotational speeds, the Simple-II model yielded more realistic 
results than the original Simple analysis at such speeds. 

Nomenclature 

General 

A gas cross-section area 

A wg gas wetted area (m 2 ) 

c average speed of molecules (m s _1 ) 

c p specific heat at constant pressure (kj kg -1 I< -1 ) 

c v specific heat at constant volume (kj kg -1 I< -1 ) 

d hydraulic diameter (m) 

Dd diameter of displacer (m) 

f rotation frequency of engine (Hz) 

/ re Reynolds friction factor 

J annular gap between piston and cylinder (m) 

I< gas conductivity 

I< g specific heat ratio [C p /C v ] 

L length of seal (m) 

Ld length of displacer (m) 

N r rotation speed of engine (RPM) 

NTU number of transfer unit 

P pressure (Pa) 

Pr Prandtl number 

Re Reynolds number 

S displacer stroke (m) 

St Stanton number 

u velocity of gas flow (m s _1 ) 

Up Piston linear velocity (m s _1 ) 

V volume (m 3 ) 

R C ond conduction resistance (m I< W -1 ) 

Greek 

p density (kg nrT 3 ) 

e regenerator effectiveness 

p, viscosity (Pa s) 


Subscript 


c 

cooler 

ck 

cooler-compression space interface 

e 

expansion space 

h 

heater 

he 

heater-expansion space interface 

hi 

hot gas input 

ho 

hot gas output 

k 

Kooler (cooler) 

la- 

Kooler (cooler)-regenerator interface 

leak 

leakage 

r 

regenerator 

rh 

regenerator-heater interface 

shuttle 

shuttle effect 

wh 

heater wall 

wk 

cooler wall 
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